Thermal concurrence mixing in a ID Ising model 
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We investigate the entanglement arising naturally in a ID Ising chain in a magnetic field in an 
arbitrary direction. We find that for different temperatures, different orientations of the magnetic 
field give maximum entanglement. In the high temperature limit, this optimal orientation corre- 
sponds to the magnetic field being perpendicular to the Ising orientation (z direction). In the low 
temperature limit, we find that varying the angle of the magnetic field very slightly from the z 
direction leads to a rapid rise in entanglement. We also find that the orientation of the magnetic 
field for maximum entanglement varies with the field amplitude. Furthermore, we have derived a 
simple rule for the mixing of concurrences (a measure of entanglement) due to mixing of pure states 
satisfying certain conditions. 

PACS numbers: 03.67.*, 75.10.Jm 



I. INTRODUCTION 

Entanglement, the rather counterintuitive non-local 
correlations exhibited by quantum systems, has recently 
become one of the most valuable resources in quantum 
information processing Q| . Over the past few years it has 
developed into a quantifiable physical resource ||, [|, ||] 
in an analogous manner to energy. Thus, the amount of 
entanglement present naturally in complex physical sys- 
tems (systems with many interacting components) now 
becomes a relevant question to ask. Condensed mat- 
ter physicists have long investigated correlations between 
parts of composite systems. Entanglement is the "quan- 
tum" or "nonlocal" part of these correlations. As such, 
it can behave very differently from the total correlations. 
For example, while correlations are averaged on mixing 
states, entanglement generally decreases. The topic of 
variation of entanglement in condensed matter systems 
with respect to the variation of external parameters such 
as temperature, field components etc., is a relatively 
unexplored and potentially rich area of study. In this 
context, as a simple initial model, Arnesen, Bose and 
Vedral have studied the variation of entanglement with 
temperature and magnetic field in a ID isotropic finite 
Heisenberg chain ||. Prior to that, Nielsen had investi- 
gated the entanglement between two qubits interacting 
via the Heisenberg interaction at a nonzero temperature 
H and O'Connor and Wootters have investigated the en- 
tanglement in the ground state of an antiferromagnetic 
isotropic Heisenberg ring Qj. In Ref. J5j, the entangle- 
ment at a nonzero temperature, being that of a thermal 
state, was called thermal entanglement. Subsequently, 
Wang has studied the quantum Heisenberg XY-model |j| 
and the two-qubit anisotropic XXZ-model j9| in a similar 
context. 

In this paper, we are going to study the thermal entan- 
glement in the ID Ising model in an external magnetic 
field j[o|. Ising- like interactions form the basic coupling 
in many proposals for experimental systems that can be 
used to perform quantum computation, see for example, 



0§§§. The ID Ising model describes a set of 
linearly arranged spins (qubits), each interacting with its 
nearest neighbors by a coupling which is proportional to 
o z ® <j z . This coupling can be diagonalized in a basis of 
disentangled states. Naively one might think that this 
implies a complete absence of entanglement in the Ising 
model. However, an external magnetic field with a com- 
ponent, however small, along a direction perpendicular to 
the z axis is sufficient to make the eigenstates entangled. 

We start in section II by considering analytically the 
case of two qubits interacting via the Ising interaction in 
a magnetic field orthogonal to the z direction. In section 



[II, we formulate a theorem for the concurrence mixing 
due to occupation of both the ground and the excited 
states. Next, in section |Tv| , we consider numerically the 
variation of entanglement with the orientation of mag- 
netic field. Before concluding, in section [y] we show that 
the kinds of behavior found for two qubits also hold for 
many qubits. 



II. TWO QUBITS WITH ISING INTERACTION 

The Hamiltonian for an isotropic quantum Ising model 
with nearest neighbor couplings in an external magnetic 
field can in the most general form be expressed as 



H = J 



<i,j> 



B 



(1) 



where the indices i, j and k label the N spins. Here we 
will consider systems in one spatial dimension with pe- 
riodic boundary conditions so the iVth spin also couples 
to the first spin. Thus we have a qubit "ring" . First we 
will consider the case of N = 2. Our Hamiltonian can 
then be written as 



H = 2 Ja z ® a z + B ■ (a ® I + 1 ® a). 



(2) 



The usual form of the Ising model has a magnetic field 
only along the z axis. This case has no entanglement 
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at all, since the Hamiltonian is diagonal in the standard 
disentangled basis {|00), |01), 1 10) , |11)}, where |0) stands 
for spin up and |1) stands for spin down. However, in this 
paper we will consider the cases when the magnetic field 
is not parallel to the z axis. 

A. Orthogonal fields 

Let us first study the special case when the magnetic 
field is perpendicular to the z axis, say B = B x. Our 
system is now described by the following Hamiltonian: 

H = 2Ja z ® a z +B (a x + a x ). (3) 

We are going to investigate the entanglement in this two- 
qubit Ising ring. In this paper, we will use the squared 
concurrence j|, [Tj|, called the tangle r, as a measure of 
entanglement. To calculate this, first we need to define 
the product matrix R of the density matrix, p, and its 
spin-flipped matrix, p — (<j y ® a v ) p* (a y <g) a y ) . Hence, we 
have 
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Magnetic field amplitude (B) 

FIG. 1: Energy levels with corresponding eigenstates, a,b,c 
and d are functions of B, and the states |$±) = (1 00) ± 
|ll))/v / 2and |#±) = (|01)±|10})/V^ are the four Bell states. 



R= P P = p{(Ty®dy)p*{ay®&y). (4) 

Now concurrence is defined by 

C = max{\i — A2 — A3 — A4, 0}, (5) 

where the A^ are the square roots of the eigenvalues of 
R, in decreasing order. In this method the standard ba- 
sis, {|00), |01), 1 10), |11>}, must be used. As usual for 
entanglement measures, the tangle ranges from (no en- 
tanglement) to 1, when the two qubits are maximally 
entangled. 

For finite temperatures we need the density matrix, 
p, for a system which is at thermal equilibrium. This 
is given by p = e~ H ' T /Z, where Z — tr(e~ H / T ) is 
the partition function (using units where the Boltzmann 
constant, ks = 1). We then solve the time inde- 
pendent Schrodingcr equation for our qubits. The en- 
ergy levels of our Hamiltonian (|^) are, in rising order, 
-2V.P + B 2 , —2 J, 2 J, 2V.P + B 2 , as in Fig. 

For zero temperature only the lowest energy level is 
populated. The tangle of this pure state can easily be 
calculated from the density matrix, for B > 0, 

j2 i 

T= j^T^ = 77^r (6) 

In Fig. H a contour plot of the tangle t, as a function of 
magnetic field amplitude, i?, and the temperature, T is 
shown. 

From Eq. (|^), it is clear that the entanglement is high- 
est for nearly vanishing magnetic fields and decreases 
with increasing field amplitude. However, Eq. (](]) is not 
valid for strictly B = 0, in which limit it seems to pre- 
dict maximal entanglement. At precisely B = 0, in fact, 
no entanglement is present (the eigenstates are same as 



Tangle t=t(T/J,B/J) 
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FIG. 2: Contour plot of the tangle of two qubits obeying an 
Ising Hamiltonian with coupling J, in a perpendicular mag- 
netic field B, for temperatures T. 



those of the usual Ising Hamiltonian without any mag- 
netic field). Hence there is a quantum phase transition 
|l0f at the point B=0 when the entanglement jumps from 
zero to maximal even for an infinitesimal increase of B. 
As we will see later this point is only one point on a 
transition line for £>-fields along the z axis. 

Let us now turn our attention to the more realistic 
case of non-zero temperatures. For a general pure state 
only one of the eigenvalues of Eq. (|4|) is non-zero and 
therefore equal to the tangle. This statement is proved 
in lemma 1 in the next section. For low temperature and 
magnetic field, i.e. B, T <C J, it is a good approximation 
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to assume that only the two lowest energy levels are pop- 
ulated. This becomes clear when we look at Fig. [| in the 
regime B« J. The lowest two levels are much closer to 
each other compared with their separation from the third 
lowest energy level (i.e. the second excited state). Thus 
when the temperature is low, only the lowest two levels 
appear in the state of the system. We will find (theorem 
1, next section) that, in our case, the combination of the 
two lowest states also combines their concurrences in the 
following weighted subtraction: 



C — max{|woC'o — w\Ci\,0}, 



(7) 



where the index refers to the ground state, while 1 refers 
to the excited state and wo and w\ are the weights of the 
ground and excited states respectively. The weights can 
be any weights from the statistics, for example Maxwell- 
Boltzmann statistics or Fermi-Dirac statistics. We call 
this concurrence mixing. In our case, the first excited 
state is the Bell state, |*") = (|01) - |10))/\/2, which 
has r = 1, and Eq. (0) reduces to 



C=\w 



J 



VJ 2 + B 2 



Wi\ 



(8) 



In general, the first term in the above equation is larger 
than the second, and in this case the concurrence de- 
creases with temperature as wq decreases and w\ in- 
creases (c/Fig. ||). In Fig. ^ we also see that, for a 
given temperature, the entanglement can be increased 
by adjusting the magnetic field and is generally largest 
for some intermediate value of the magnetic field. This 
effect can be understood by noting that wq increases with 
increasing B as the energy separation between the lev- 
els increase, but J/V J 2 + B 2 decreases. As a result the 
combined function reaches a peak as we vary B and de- 
creases subsequently, inducing analogous behavior for the 
concurrence. 



III. CONCURRENCE MIXING 

In this section, we are going to formulate and prove 
a useful concurrence mixing theorem. We begin with a 
lemma which illustrates the method used in the theorem. 
The results of lemma 1 appear in Ref . || . 

Lemma 1: Let p be a pure density matrix. Then the 
product matrix R = pp, where p = (a y ® &y)p*(cr y £g> <j y ), 
has only one non-zero eigenvalue, and its value is the 
concurrence squared, i.e. the tangle. For a general pure 
state, \a) = a|00) +6|01) + c| 10) + rf|ll) the concurrence 
is C = 2\ad— bc\ or, written as a Schmidt decomposition, 
C = 2\cqCi\, where Ci are the two Schmidt coefficients. 

Proof: Consider a general pure density matrix p = \ot)(a\. 
By writing out the product matrix 



we see directly that |mo) = \a) is an eigenstate with 
eigenvalue Kalcr,, <g> <7j,|a)| 2 . It is always possible to 
find three more vectors \a^),k — 1,2,3 all of them 
linearly independent of each other and orthogonal to \a). 
Thus the remaining three eigenvectors can be written 
as | Ufe) = cT y ® d-y\a^),k — 1,2,3, all with eigenvalue 



zero. From Eq. (|5|) we now get the concurrence to be 
C = \(ot\a y ® Vy\a)\ = 2\ad — bc\. Thus, for a pure state, 
concurrence can be defined as the absolute expectation 
value of the operator a y <g> a y . By tracing out one qubit 
and solving for the eigenvalues, which are equal to the 
square of the Schmidt coefficients, c;, of the remaining 
density matrix, we find that the concurrence also can be 
written as C = 2\cqCi\. 

Theorem 1: Consider two pure states of the same 
system \a m ) and \a n ). If the spin-flip overlap is zero, 
i.e. 



(a m \o-y ® cr y \a n } 



0, 



(10) 



then the concurrence of the mixture of the two pure states, 
with weights Wi, can be expressed as 



C mixed — l^mCm ^nCn|- 



(11) 



Proof: Let pi = \a.i)(a.i\, i = m,n, be our two 
pure states. From our lemma, we have 

Pip~i\u i0 ) =C?\uiq) (12) 
PiPi\Uik) =0, A; = 1,2,3 (13) 

where \uio) = \oti) and \uik) = <r y ® <J y \a^ k ). Let us write 
our mixed state, p, as a weighted average of the pure 
density matrices, p = w m p m + w n p n - Since p is only a 
linear transformation of p, we also have p = w m p~ m + 
w n p~ n . Using these assumptions, we can write down the 
product matrix 

PP = W^Pmp'm + W m W n (p m p~ n + p n P~m) + W 2 n p n p n . 

(14) 

Our condition ( jlpj ) makes the cross terms drop out, since 
PiPj = \ai){ai\a y ® a v \aj){aj\a y ® a y = 0. (15) 



pp = \a)(a\a y <g> a y \a)(a\a y (g) a y , 



(9) 



Furthermore, condition (jlfj) together with Eq. (|12j) gives 
the following relation 

PiPi\ujo) = \ai){ai\a y <® a v \cti)(ai\(j y ® a v \ctj) 

= S lJ C 2 \u j0 ). (16) 

The Eqs. (|l4|)-(]l6|) give two of the four eigenequations of 
the product matrix 

pp\vi) = (5 im w 2 m Cl l + 5 in wlCl)\vi) , (17) 

where \vi) = \uio) — \on). Since these two eigenvectors 
only span two dimensions in the four-dimensional space, 



4 



we can always find another two vectors \a^),k = 2,3 
which are linear independent of each other and or- 
thogonal to the two eigenstates. Thus the last two 
eigenvectors are \vf.) = fry ® a y \a^),k = 2,3, both 
with zero eigenvalue. Eq. (^) now gives our mixed 
concurrence formula. 

Theorem 1 applies to any system with a mixture 
of two pure states satisfying condition ( JTo| ) . It can easily 
be extended to apply to the mixing of more pure states. 
The requirement is then that condition ( p0| ) must hold 
for all pairs of pure states. 

In our case, it would have been interesting if, when 
including all four levels, the concurrence could be calcu- 
lated as 

C = max{2w k C k - w i C i> °>- ( 18 ) 

i 

In fact, for a three- level approximation (involving the 
first three levels), this formula is correct, however the 
exact four-level concurrence is not in agreement with Eq. 
( |l8| ) , because condition ([Io|) does not hold for mixing the 
ground state with the third excited state. 

IV. ARBITRARY FIELDS 

In section II we treated a case of a magnetic field or- 
thogonal to the Ising direction. We are now going to 
generalize this to arbitrary magnetic fields. The new 
Hamiltonian can be written 

H = 2J & z ® o z 

+Bsin9 (o x ® I + I <g> cr x ) 
+Bcos6(a z ®i + i®& z ), (19) 

where 9 is the angle between the magnetic field and the 
Ising direction. It is sufficient to consider variation of B 
in a plane containing the Ising direction, because in 3 
spatial dimensions the Hamiltonian possesses rotational 
symmetry about the z axis. 

The expression for the tangle is analytically solvable. 
However, because of a difficult cubic equation in the 
diagonalization, the expressions are complicated, so we 
present the results in graphical form. At zero tempera- 
ture, Fig. shows our solution when the tangle is plotted 
as a function of B x — B sin 9 and B z = B cos 9, and Fig. 
H shows the solution when the tangle is plotted as a func- 
tion of the amplitude B and angle 9. 

We notice that the region around B x = for all 
\B Z \ < 2 J has the highest possible entanglement. At 
exactly at B x = there should not be any entangle- 
ment, the white region of Fig. || indicates a quantum 
phase transition at B x = (the sharpness of the tran- 
sition being illustrated by the fact that the zero entan- 
glement line at B x — is so thin that it is invisible). 
For small angles, 9, there are two energy levels close to 
the energy — 2J with corresponding states close to the 



Tangle t=t(B /J,B /J) 




FIG. 3: Contour plot of the tangle at zero temperature in 
a Cartesian coordinate system. Note that the line of non- 
entangled states at B x = for \B Z \ < 2 J can not be seen. 
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FIG. 4: Contour plot of the tangle at zero temperature in a 
spherical coordinate system. 

Bell states |* ± ) = (|01) ± \10))/y/(2). Thus we get a 
maximally entangled qubit pair in the limit. However, 
at B = B z z, with \B Z \ < 2 J, the states are degenerate 
with no entanglement as a result. In the case B z > 2,1 
(B z < —2 J) the ground state is always the non-entangled 
state 1 11) (| 00)). Also notice that 9 = it/2 corresponds 
to our earlier orthogonal case, thus the tangle follows Eq. 
(^). Even when \B X \ is increased to the point where it 
starts to dominate, the spins will simply align along B x 
and give a disentangled state. Thus the entanglement 
falls off with increasing strength of the magnetic field in 
either direction. 

Let us now look at the case of finite temperature (ther- 
mal entanglement). The first excited state is l^ - ) which 
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FIG. 5: Contour plots of the tangle in Cartesian coordinates 
for various finite temperatures, T. 



lies at the energy —2 J. This state is totally independent 
of magnetic field, thus the tangle corresponding to this 
state forms a constant plane at 1 . Fortunately, condition 
( |iTj| ) in our theorem is also satisfied, which makes Eq. (Q) 
valid. In Figs. || and|| a numerical solution is shown at a 
low temperature. Note how fast the tangle drops to zero 
for a low |B x |-component. This does not contradict the 
concurrence mixing formula, as the weights Wi also de- 
pend on the magnetic field through the energy. The fast 
drop in the tangle is due to the degeneracy at B x = 0. 
The smallness of the energy difference at low values of 
\B X \ makes the two levels almost equally populated even 
for small temperatures. The line of zero entanglement 
at B x = for T = has broadened into a region of 
almost zero entanglement in the finite temperature case 
(compare Figs. ||, ||). 

It is also interesting to see that there exists an angle 
9 = 6*(B, T), where the entanglement is maximum for a 
given temperature and amplitude 



maxr(B,T) = t(9*(B,T)). 



(20) 



This feature can be explained heuristically if we assume 
that with B x and B z fixed, the entanglement should 
change continuously with temperature. We know that 
increasing the temperature widens the low entanglement 
zone around B x = and the entanglement has to fall off 
for large \B X \. So it is expected that at some intermedi- 
ate value of \B X | (and hence 9) the maximal entanglement 
will be reached. As we increase the temperature further, 
the near-zero entanglement zone centered around B x = 
widens even more and pushes the entanglement maxima 
away to higher and higher values of \B X \. The highest 
value of the tangle tends more and more towards orthog- 
onal fields (c.f. Figs. |[ ||). The preferred angle traverses 
from 9 = at zero temperature (Fig. ^) to 9 — n/2 at 
T w J (Fig. ^J). In the classical limit, i.e. at very high 
temperatures all entanglement fades away as expected, 
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FIG. 6: Contour plots of the tangle in spherical coordinates. 
The angle, resulting in the maximum entanglement for a given 
magnetic field amplitude, varies as a function of the amplitude 
and temperature. 



because the state is completely mixed. 

We can use the two-level approximation to get an es- 
timate of the angle creating maximal entanglement. Let 
us first estimate the ground state energy, e. Since we 
know that the lowest two energy levels are very close and 
the first excited state is —2 J, we can use the approxima- 
tion —2 J + e « — 4J while solving for eigenvalues of the 
Hamiltonian (|l9|). The ground state energy is then 



-2 J 



ABi 



4 J 2 - Bl 



J. 



(21) 



For the above approximation to hold, we must operate 
in a region not too close to the poles at B z — ±2J. 
This gives the following energy difference between the 
two lowest levels 



Ae = -2J-e = 



ABl 



AJ 2 - B1 



J. 



(22) 



The two terms of the ground state energy (21) can be 
considered as the first two terms in a Taylor expansion. 
This suggests that we instead write 



'=-^ /1 + T = -W 1 + 47^' 



(23) 



Let us, from here onwards, measure all energies in units 
of J, i.e. let J = 1. For a magnetic field only in the x 
direction, the concurrence of the ground state is Cq = 
t 1 / 2 = — 2J/e, where r is given by Eq. For non-zero 
B z this formula remains an excellent approximation (as 
we have verified numerically). Substituting Eq. ( p3[ ) for 
e gives 



Cn 



1 



VT+Ae 



(24) 
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Recall that the first excited state is \^^) with concur- 
rence C\ = 1. From the concurrence mixing theorem Jll| ) 
we get an approximation of the thermal concurrence. If 
Wi are weights following Maxwell-Boltzmann statistics, 
then the maximum entanglement is reached when the 
following condition holds: 



T{1 + e Ae / T ) = 2(1 + Ae + vT+"A7). 



(25) 



In the region we are interested in, where temperature is 
not bigger than ~ 20% of the coupling constant, J, the 
temperature is also much smaller than Ae, and Eq. (I2F 
can be further simplified to 



Ae = Tin- 



T 



(26) 



Remember that the ground state energy is still a function 
of the magnetic field. In order to fix the temperature, Ae 
has to be fixed, i.e. maximum entanglement is reached 
in the cross section between the energy surface and the 
constant energy plane that follows from Eq. (p6f). From 
Eq. (^2j) it is clear that this is described by an ellipse in a 
B x B z -plane. Using the field amplitude B as a parameter, 
the angle we defined in Eq. (E0) is given by 



sm( 



= ± 



Ac 



B 2 



4 - Ae B 2 



(27) 



assuming that the parameter \B\ > ^/e. 

Another way to parameterize the optimum line is to 
let B z be the parameter and solve for B x , 



B x =±1 /At 



B 2 
1 z - 

4 



(28) 



with \B Z \ < 2. Again keep in mind that Eqs. fl27j ) and 
( p8| ) follow from the assumption that Ae « 0, and there- 
fore the ellipse is not closed around the poles at B z = ±2 
(cf Figs. |0). 



V. QUBIT RINGS 

As the expressions for the thermal tangle of the gen- 
eral two-qubit case are already quite complicated, we can 
not expect to find any easily manageable analytic expres- 
sions in the many-qubit case. Instead, we have performed 
numerical simulations, which gives the entanglement be- 
tween neighboring qubits as shown in Fig. [7]. 

The behavior for even N rings is quite similar to that 
of the two qubit case. To understand why there is an 
extra low entanglement zone around B z = in the case 
of odd N rings, one has to go back the basic cause for 
entanglement arising in the Ising chain. It results from 
the competition between the term a z <g> a z + B z a z trying 
to impose spin order in the z direction and B x a x try- 
ing to impose spin order in the x direction. In the odd 
qubit case, it is impossible for all neighboring spins to be 



Tangle x, N = 3, a = 1 
2| 



Tangle x, N = 4, a = 1 




Tangier, N = 5, a = 1 
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FIG. 7: Contour plots of the tangle at temperature T = 
0.10 J. N is the number of qubits in the chain and a=l 
means that entanglement is measured between two neighbor- 
ing qubits. 



oriented oppositely, so the ordering power of a z ® a z is 
significantly reduced. In such circumstances, it is mainly 
the competition between B z a z and B x a x (albeit aided 
by the small a z ® a z interaction) which determines the 
entanglement. Thus the high entanglement values near 
B z — present for the two qubit (and all other even N) 
cases disappear for odd N. Note also that the entangle- 
ment in the odd N case is somewhat larger in magnitude 
compared with the even N case. This is a result of the 
fact that the two terms a z ®a z and B z a z compete for the 
type of z ordering (parallel or anti-parallel neighboring 
spins). In even N case, this competition is much stronger 
and this tends to lower the entanglement by reducing the 
net effect of z ordering terms with respect to x ordering 
terms. As the number N of qubits in the chain is in- 
creased, the difference between even and odd N chains 
should disappear (because for large N, adding or remov- 
ing an extra qubit from the chain should not make a 
significant difference). This effect is clearly seen in Fig. 

where the difference in appearance between the plots 
for N — 3 and N — 4 is much greater than that between 
N = 5 and N = 6. 







Entanglement can also be calculated between non- 
neighboring qubits with the results shown in Fig. || for 
next-nearest neighbors. Again we observe that the even 
N case has lower entanglement on average than the odd 
N case. Also in the odd qubit case, the entanglement 
between next-nearest neighbors is somewhat complemen- 
tary to that between nearest neighbors (this can be seen 
for example by placing the plots for N = 5 in the two 
cases on top of each other). Thus the amount of entan- 
glement between pairs of nearest neighbors and pairs of 
next-nearest neighbors can be controlled by varying the 
field direction. 
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Tangle x, N = 4, a = 2 Tangle x, N = 5, a = 2 




Tangier, N = 6, a = 2 Tangle x, N = 7, a = 2 




FIG. 8: Contour plots of the tangle at temperature T = 
0.10 J. N is the number of qubits in the chain and a=2 means 
that entanglement is measured in pairs with one qubit in- 
bctwccn. 



qubit analytically and three through seven qubits numer- 
ically. One of the most interesting results is the fact that 
for a given temperature, the (nearest neighbor, pairwise) 
entanglement in the ring can be maximized by rotating 
the magnetic field (at fixed magnitude) to an optimal di- 
rection. This can be regarded as magnetically induced 
entanglement. The pairwise entanglement between next- 
nearest neighbors can be maximized by rotating the field 
to a different direction. We have also proved a theorem of 
mixing of concurrences which is applicable to any system 
in which the pure states in the mixture have no spin-flip 
overlap. 

So far we have have only considered pairwise entangle- 
ment. In future work we will estimate the entanglement 
between three or more qubits in the ring and also focus 
on investigating ways to detect the natural entanglement 
in Ising models, and on investigation of the entanglement 
in the large variety of available condensed matter models 
of interacting systems. 



VI. CONCLUSIONS 
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